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Magnetic properties of the S — 1/2 antiferromagnet (a-Cu 2 V 207 have been studied using magneti¬ 
zation, Quantum Monte Carlo (QMC) simulations, and neutron diffraction. Magnetic susceptibility 
shows a broad peak at ~ 50 K followed by an abrupt increase indicative of a phase transition to 
a magnetically ordered state at Tn = 33.4(1) K. Above Tat, a fit to the Curie-Weiss law gives a 
Curie-Weiss temperature of © = —73(1) K suggesting the dominant antiferromagnetic coupling. 

The result of the QMC calculations on the helical-honeycomb spin network with two antiferromag¬ 
netic exchange interactions Ji and J 2 provides a better fit to the susceptibility than the previ¬ 
ously proposed spin-chain model. Two sets of the coupling parameters Ji : J 2 = 1 : 0.45 with 
Ji = 5.79(1) meV and 0.65 : 1 with J 2 = 6.31(1) meV yield equally good fits down to ~ Tn- 
Below Tn, weak ferromagnetism due to spin canting is observed. The canting is caused by the 
Dzyaloshinskii-Moriya interaction with an estimated 5c-plane component \Dp\ ~ 0.14Ji. Neutron 
diffraction reveals that the S = 1/2 Cu^^ spins antiferromagnetically align in the Fd'd'2 magnetic 
space group. The ordered moment of 0.93(9) /is is predominantly along the crystallographic a-axis. 

PACS numbers: 75.30.Cr, 75.25.-j, 75.10.Pq, 75.30.Cr, 75.50.Ee 


I. INTRODUCTION 

Low-dimensional quantum magnetism has attracted 
much interest from both theoretical and experimen¬ 
tal condensed matter physicists for many decades. 

It is known that in a one dimensional {ID) antiferro¬ 
magnetic system, long range order is absent even at 
zero temperature,^ leading to various fascinating mag¬ 
netic ground states and phenomena at low tempera¬ 
tures such as the spin-Peierls state in GuGeOs^’^ and 
TiOGl,^’^ the singlet ground state in the alternating spin- 
chain system (VO) 2 P 207 ,^ the fractional spinon excita¬ 
tions in GuS 04 ' 5 D 20 ,^ and the Bose-Einstein conden¬ 
sation (BEG) of magnons in the double spin-chain sys¬ 
tem TlGuGls,^^ in which weakly interacting dimers are 
formed at low temperatures and BEG is realized as the 
field-induced 319 magnetic ordering. However, some ap¬ 
proximate ID antiferromagnets exhibit long-range order 
with a remnant of quantum fluctuations in a form of 
quantum renormalization of spin waves.Hence, in or¬ 
der to apprehend diverse physics of these low-dimensional 
quantum magnets, it is crucial to identify a spin network 
and relevant underlying interactions that consequently 
cause magnetic ordering and govern spin dynamics. 

Among various compounds, copper-based oxides with 
Gu^+ ions {3d^), titanium-based oxides with Ti^+ 
ions (3(1^), and vanadium-based oxides with ions 
{3d^) are generally considered as a good realization of 
the low-dimensional spin- 1/2 system. Gopper divana¬ 


date GU 2 V 2 O 7 is a promising realization of the low¬ 
dimensional spin-1/2 system. There are two polymorphs 
of GU 2 V 2 O 7 that are generally found in a stable phase, 
namely the a-phase and the /-phase with a structural 
phase transition at around 605°G.^^ Other forms such 
as the 7 - and /'-phases were also reported as complex 
and unstable high-temperature phases.The a-phase 
of GU 2 V 2 O 7 crystallizes in the orthorhombic system of 
space group Fdd2 with a = 20.645(2) A, 5 = 8.383(7) A, 
and c = 6.442(1) Structurally, Gu^+ ions that are 

surrounded by five oxygen ions appear to form a chain 
of edge sharing GuOs polyhedra while the nonmagnetic 
ions, each of which is surrounded by four oxygen 
ions to form V 2 O 7 double corner-sharing tetrahedra, are 
intercalated between the chains as shown in Eig. 1(a). 
Within the 6 c-plane, the S' = 1/2 spins of Gu^+ ions form 
a zigzag chain and interact with their nearest neighbors 
via two inequivalent Gu-O-Gu pathways. The distance 
between the two copper ions Gul and Gu4 connected 
by the blue bonds in Eig. 1 (a) is approximately 3.1 A. 
An intra-chain interactions between these two coppers 
are formed by the bridging of Gul(Gu 2 ) with Gu4(Gu3) 
along [oil]([oil]) via 0(2) and 0(3). The deviation of 
the Gu-0(2)-Gu angle from 90° (Eig. 1 (b)) leads to a 
preferable and strong antiferromagnetic interaction along 
the chain. In the previous study, this system has been 
proposed as a realization of a zigzag spin-chain model. 
However, the possibly strong next-nearest-neighbor in¬ 
teractions between the Gu^+ spins, which link the zigzag 
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FIG. 1. (Color online) (a) Crystal structure of a-Cu 2 V 207 . 
The chains of edge sharing CuOs polyhedra are along [Oil] 
and [Oil] with two corner-sharing V 2 O 7 tetrahedra interca¬ 
lated between the chains. Blue (Ji) and green (J 2 ) bonds rep¬ 
resent the antiferromagnetic couplings between the nearest- 
neighbor, intra-chain spins and next-nearest-neighbor, inter¬ 
chain spins, respectively, (b) The DM interactions, the direc¬ 
tions of which can be determined from the cross product of 
X Yij and x Vij for two inequivalent Cul-0(2)-Cu4 and 
Cul-0(3)-Cu4 bonds, respectively, yield a non-zero value, 
(c) On the other hand, the equivalent Cul-0(1)-Cu2 path¬ 
ways connecting the inter-chain spins give rise to a net van¬ 
ishing DM interaction. The structures are visualized using 
VESTA}^ 


chains along the crystallographic a-axis as shown by the 
green bonds in Fig. 1(a), suggest a highly anisotropic 
and non-trivial spin network with three coordinate spins 
for each site that resembles a disconnected honeycomb 
(helical-honeycomb) lattice when viewed along the 6 -axis 
(Fig. 9). The next-nearest-neighbor bond of length 4.0 A, 
which is slightly longer than the nearest-neighbor bond, 
is formed by the bridging of Cul(Cu3) with Cu2(Cu4) 
via two equivalent 0(1) ions (Fig. 1 (c)) with the Cu- 
0(1)-Cu angle of 104°, possibly leading to an antifer¬ 
romagnetic interaction that could be comparable to the 
intra-chain interaction. This sizable inter-chain interac¬ 
tion that enables the forming of the helical-honeycomb 
lattice is not surprising. Using first-principle calculations 
on its cousin phase / 3 -CU 2 V 2 O 7 , Tsirlin et also sug¬ 
gested that the ;5-phase could be better described by a 
honeycomb lattice than by a spin-chain model. The co¬ 
ordination number of three is the lowest for any 2D lat¬ 
tices and only d. ID chain has a lower coordination num¬ 
ber. Therefore, the helical-honeycomb model could share 
many properties unique to the low-dimensional magnets 
even though the lattice extends in three dimensions. 

Magnetic susceptibility and heat capacity on a pow¬ 
der sample reveal that the antiferromagnet (a-Cu 2 V 207 
magnetically orders at a Neel temperature T/v of 
34 The magnetic ground state of (a-Cu 2 V 207 

was previously proposed to be a canted antiferromag¬ 
netic spin-chain.The canting of spins is due to 


Dzyaloshinskii-Moriya (DM) interations, which give rise 
to the presence of weak ferromagnetism below Tyr. How¬ 
ever, the magnetic structure and spin dynamics of a- 
CU 2 V 2 O 7 have not been studied. Among the various 
probing techniques, neutron scattering is the most pow¬ 
erful at revealing the microscopic properties of such mag¬ 
netic materials; however, a good quality large single crys¬ 
tal is required. Here we report the first detailed study of 
magnetic properties on single-crystal a-Cu 2 V 207 using 
magnetization. Quantum Monte Carlo (QMC) simula¬ 
tions, and neutron diffraction measurements. 

The paper is organized as follows. In Sec. H, we 
briefly discuss single-crystal growth and experimental 
techniques used to characterize and study magnetic prop¬ 
erties of the compound. The results in Sec. HI are divided 
into three parts, where x-ray diffraction (Sec. HI A), mag¬ 
netization and QMC simulations (Sec. HIB), and neu¬ 
tron diffraction (Sec. HIC) will be discussed. Finally, we 
end with the summary in Sec. IV. 


II. EXPERIMENTAL DETAILS 

Prior to the single-crystal growth, powder CU 2 V 2 O 7 
was prepared from high purity CuO and V 2 O 5 . The 
chemicals were dehydrated and weighed with stoichio¬ 
metric ratio and then ground thoroughly with ethanol. 
The mixture was then calcined at 500°C for 24 hours. 
The obtained powder was inserted into a quartz tube. 
The bottom end of the tube was shaped into a taper for 
seed selection while the top end was tightly closed with 
silica wool. The sample was melted in air at 850° C for 10 
hours to ensure homogeneity and then lowered through a 
constant temperature gradient of about 40°C/cm inside 
a vertical Bridgman furnace at a speed of 1 cm/day. The 
sample was finally cooled from 700° C to room tempera¬ 
ture at a rate of 5°C/min. Single crystals were extracted 
from the quartz tube by mechanical separation. A pure- 
phase powder sample of a-Cu 2 V 2 07 was also synthesized 
by the standard solid state reaction and used in powder 
neutron diffraction measurements. 

Small pieces of the crushed crystals were collected and 
ground for powder x-ray diffraction measurements us¬ 
ing CviKa radiation. The results were fit using the Ri- 
etveld method in FullProf.^^ To confirm the crystal struc¬ 
ture, single-crystal x-ray diffraction data were collected 
at room temperature using a Bruker X 8 APEX CCD 
diffractometer with MoKa radiation. The refinements 
were done using the software ShelXle.^^ The magnetiza¬ 
tion M of the single-crystal sample was measured to the 
lowest temperature of 2 K using a superconducting quan¬ 
tum interference device (MPMS-XL, Quantum Design). 

QMC simulations with the LOOP algorithm^^ were per¬ 
formed using the simulation package ALPS.^^ The mag¬ 
netic susceptibility was calculated on a cluster of 100 
spins for the spin-chain model and up to 432 spins (27 
unit cells) for the helical-honeycomb model with a pe¬ 
riodic boundary condition in the temperature range of 
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0.01 < t < 5 (t = ksT/J) using 100 000 sweeps for ther- 
malization and 500 000 Monte Carlo steps after thermal- 
ization. 

To check crystallinity quality and investigate the mag¬ 
netic transition, single-crystal neutron diffraction was 
performed at the BT7 Double Focusing Thermal Triple 
Axis Spectrometer^^ at NIST Center for Neutron Re- 
seach, USA, on a single crystal with a mosaic of 0.8°. 
Elastic neutron scattering were performed at 2.5 K and 
50 K using fixed incident energies of 14.7 meV and 
30.5 meV. The position-sensitive detector (PSD) was em¬ 
ployed in a two-axis mode with open - 80' - sample - 
- 80' (radial) - PSD horizontal collimations to map out 
the broad reciprocal space in the (hkO) scattering plane 
(Fig. 6(a)). A detailed investigation of the nuclear and 
magnetic Bragg reflections were performed using a triple¬ 
axis mode with a single detector and the horizontal col- 
limations of open - 80' - sample - 80' - 120'. For all 
diffraction measurements, one pyrolytic graphite (PC) 
filter and two PC filters were placed along the incident 
and scattered beams, respectively, to reduce higher-order 
neutron contamination. Neutron diffraction on the pow¬ 
der sample were performed at the high-resolution neu¬ 
tron diffractrometer Echidna, ANSTO, Australia using 
neutrons with wavelength 2.44 A(13.7 meV). The data 
were collected at 3 K and 40 K to extract the magnetic 
scattering. Additional single-crystal neutron diffraction 
measurements were carried out at 4 K and 50 K us¬ 
ing the Laue diffractometer Koala, ANSTO, Australia. 
The Laue data images were processed using the LaueG 
software. 


III. RESULTS AND DISCUSSION 
A. X-ray diffraction 

The obtained single crystals are shown in the inset of 
Fig. 2. The largest crystal was about 1 x 1 x 0.5 cm^ 
with a mass of 1.4 g. The naturally cleaved facet is the 
( 1 , 0 , 0 ) plane. Small crystals were collected for the x- 
ray diffraction measurements while large crystals were 
allocated for the neutron scattering experiments. The 
result of the powder x-ray diffraction (Fig. 2) shows that 
the major phase of the crystals is (a-Cu 2 V 2 07 (^95%) 
with a small amount of impurities that can be identi¬ 
fied as / 3 -CU 2 V 2 O 7 and Cuo. 64 V 205 .^^ We note that no 
trace of these impurities is detected in the powder sam¬ 
ple (Fig. 7). The lattice constants obtained from the 
refinement are a = 20.678(6) A, 6 = 8.405(2) A, and c 
= 6.446(2) A, which are in good agreement with those 
reported in Ref. 18. The impurities are still present af¬ 
ter several crystal growth attempts with different cool¬ 
ing conditions, which are an important factor to con¬ 
trol the ratio of a and P phases. We found that the 
fraction of /3—phase increases with an increasing cool¬ 
ing rate. Hence it is crucial to slowly cool the sample 
through the phase transition temperature at 605° C to 
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FIG. 2. (Color online) Powder x-ray diffraction was mea¬ 
sured on crushed single crystals, black for the observed data, 
red for the calculated pattern, blue for the difference, and 
green for the Bragg positions. The symbols x and o indicate 
CU 0 . 64 V 2 O 5 and / 3 -CU 2 V 2 O 7 impurities, respectively. The in¬ 
set shows a photograph of the 0 -CU 2 V 2 O 7 single crystals. 


avoid a mixture of the two polymorphs. It should be 
noted that these impurities comprise only small percent¬ 
age, are most likely in a powder form, and hence will 
not mislead the interpretation of the neutron scattering 
data (Figs. 6 and 8). Room-temperature single-crystal 
x-ray diffraction was performed on the crystal with a 
few hundreds of micrometers in size. The data were re¬ 
fined against space group Fdd2 with the previously re¬ 
ported lattice parameters'^ yielding the agreement fac¬ 
tor Ri = 0.039 for 1031 reflections with Fobs < 4cr(Fo6s). 
The result is shown in Table 11. The refinement result 
from single-crystal neutron diffraction measured at 50 K, 
which will be discussed later, is also shown in the table 
for comparison. 

B. Magnetic susceptibility 

To investigate the magnetic transition on single-crystal 
a-Cu 2 V 207 , we measured the magnetic susceptibility 
(x = M/H) as a function of temperature when the ap¬ 
plied magnetic field of 1 T was parallel and perpendicular 
to the a-axis in the zero-field-cooled mode. In Fig. 3(a) 
the susceptibility exhibits clear anisotropy at low tem¬ 
peratures and shows a sharp transition at ^ 33 K in 
agreement with Tjy = 33.4(1) K obtained from the order 
parameter measured by neutron diffraction on a single 
crystal at BT7 (the inset). The Neel temperature is con¬ 
sistent with that obtained from the susceptibility mea¬ 
sured on our powder sample (not shown), as well as with 
those from the previous powder-sample studies. 

Above 100 K, a linear fit of as a function of T to 

the Curie-Weiss law shown in Fig. 3(b) 
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FIG. 3. (Color online) (a) Magnetic susceptibility was mea¬ 
sured with the applied field of 1 T parallel (H\\a) and per¬ 
pendicular {H±a) to the crystallographic a-axis. The inset 
shows an order parameter as a function of temperature mea¬ 
sured by neutron scattering at the magnetic Bragg reflection 
(0,2,0). (b) The Curie-Weiss law (red) is fit to the suscepti¬ 
bility with i4Ta. 

gives a Curie-Weiss constant C = 0.545(2) cm^K/mol Cu 
and a Curie-Weiss temperature 0 = —73(1) K. The neg¬ 
ative Curie-Weiss temperature suggests that the domi¬ 
nant exchange interactions are antiferromagnetic. From 
the Curie-Weiss constant, the calculated effective mo¬ 
ment jiieff = y^^ksC/NA = 2.087(4) is obtained. 
This value is slig htly high er than the spin-only value 
of jiieff = gi^B\/S{S + 1) = 1.73 iiB for ^ = 2 and 
S = 1/2 for Cu^+ ions. The order of frustration de¬ 
fined by / = |0/T/v| is ^ 2.2, which suggests that 
the spin interactions are not strongly frustrated (a typ¬ 
ical value for strongly frustrated systems is / > 10). 
Assuming the mean field approximation, one can cal¬ 
culate the antiferromagnetic exchange interaction from 
0 = —zJcwS{S + 1)/3/cb, where z is a number of the 
coordinate spins and S = 1/2. The calculations give 
Jew = 12.6(2) meV for the spin-chain model with z = 2 
and 8.4(2) meV for the helical-honeycomb lattice, in 
which each spin has three coordinate spins, two along 
the zigzag chain and one between the chains (Figs. 1(a) 
and 9(a)), giving z = 3. 

Below T/v, when the magnetic field is applied parallel 
to the a-axis, a small cusp can be observed at the mag¬ 
netic ordering transition. This cusp is a signature of an 
antiferromagnetic transition and suggests that the spins 
align anti-parallel along the crystallographic a-axis. On 


the other hand, when the field is applied perpendicular 
to the a-axis, the susceptibility shows a broad maximum 
around 50 K suggesting a rise of short-range correlations, 
typical for low-dimensional magnets, before an abrupt in¬ 
crease indicative of long-range ordering at lower temper¬ 
atures. The weak ferromagnetism for T a below T/v is 
a result of small spin canting due to the DM interactions. 
Hence, to first approximation, the spin Hamiltonian can 
be described by 

U = J2 • S, + D,,. • (Si X S,)], (1) 

(bi> 

where Jij denotes the exchange couplings representing 
the nearest-neighbor, intra-chain Ji and next-nearest- 
neighbor, inter-chain J 2 interactions as shown by the blue 
and green bonds in Fig. 1(a), respectively. Dij is the DM 
vector whose strength is proportional to spin-orbit cou¬ 
pling and scaled with the exchange interaction between 
Si and Sj. The DM interaction is present in (a-Cu 2 V 207 
for the Cul-Cu4 bond (Fig. 1(b)) since there is no in¬ 
version center between the magnetic ions.^^ The 

magnitude of spin canting 77 resulted from the DM in¬ 
teractions was estimated to be 2° from the previous 
magnetization measurements on the powder sample. 

To determine the spin canting and DM parameter on 
the single-crystal sample, the magnetization measure¬ 
ments as a function of magnetic field up to a maximal 
field of 5 T were performed with two orthogonal mag¬ 
netic field orientations, i.e. H || a and H ± a. Figure 4 
shows the magnetization measurement on the single crys¬ 
tal with the applied magnetic field perpendicular to the 
a-axis. Above T/v (T = 50 K) the magnetization is linear 
throughout the measuring field range. However, below 
T/v (T = 2 K) starting from zero applied field the mag¬ 
netization sharply rises to a finite value with only a slight 
increase in magnetic field before attaining the same lin¬ 
ear response as that observed in the 50 K data. The 
rapid increase of the magnetization, which is indicative 
of weak ferromagnetism due to the spin canting, is not 
observed for || a as shown in the lower inset. The small 
kink around zero field is most likely due to slight mis¬ 
alignment. The same measurements on the power sample 
show a similar weak ferromagnetic component as shown 
in the upper inset. However, the jump is less sharp and 
about a factor of two smaller due to powder average. For 
the single-crystal data, a hysteresis loop, which is typical 
for ferromagnetism, is not clear with a very small coercive 
field (< 10 Oe, which is the resolution of the measure¬ 
ments), but it is more pronounced for the powder data, 
possibly due to the powder average over all orientations, 
which could broaden the magnetization reversal. 

Quantitatively, the red line in Fig. 4 (and in the up¬ 
per inset for the powder sample) denotes a linear fit to 
the magnetization for H > 5000 Oe and is extrapolated 
to intercept the 7 /-axis to obtain the value of M(0), the 
canted moment at zero field. From the single-crystal 
(powder) data, M(0) = 0.0698(1) i^b (0.0364(1) jib) is 
obtained. On the other hand, with H \\ a the value of 
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FIG. 4. (Color online) Magnetization as a function of mag¬ 
netic field was measured on single-crystal a-Cu 2 V 207 with 
i/_LaatT = 2K and 50 K. The lower inset shows the data 
for H II a. The upper inset shows the magnetization as a 
function of magnetic field measured on a powder sample at T 
= 5 K. 


M( 0 ) is approaching zero as expected since the predomi¬ 
nant spin component is anti-parallel along the a-axis and 
the canted moments only stay within the 6 c-plane. We 
note that there exist two sets of the zigzag chains along 
[Oil] and [Oil] on the alternating planes. These chains 
are about 75° with respect to each other (Fig. 9(b)). We 
will assume that for the single-crystal data the measured 
canted moment M( 0 ) is a saturated value, where the 
canted moments on both sets of the chains are aligned 
along the direction of the applied magnetic field. The 
canting angle can be calculated from the relation r] = 

sin“^ ^ ) • Given the expected spin-only S = 1 / 2 , 

^ = 2 , and M(0) from the single-crystal data, the canted 
moment of 77 = 4.0° is obtained. 

The direction of the DM vector is determined by 
X Vij, where Vij is a unit vector connecting the spins 
Si and Sj, and is perpendicular to Vij and points to¬ 
ward the oxygen ligand as shown in Fig. 1(b). The two 
pathways Cul-0(2)-Cu4 and Cul-0(3)-Cu4 are struc¬ 
turally inequivalent giving rise to non-compensating DM 
interactions. We note that the DM interactions arising 
from the two equivalent Cul-0(1)-Cu2 bonds (Fig. 1(c)) 
compensate each other resulting in the vanishing DM vec¬ 
tor. Although the resulting x does not restrict the 
DM vector to a specified high-symmetry plane, the cross 
product of a spin pair, x S^, constrains the relevant 
component of the DM vector, which causes the canting, 
to be only within the 5c-plane. We will later show in 
Sec. IIIC that the canted moments in the 5c-plane are 
parallel, and hence the component of x Sj along the 
a-axis vanishes. If only the relevant interaction Ji be¬ 
tween Cul and Cu4 is taken into account, the in-plane 
DM parameter Dp can be related to the canting angle 77 


and Ji through the following relation: 


tan(277) = 


Dp 

Ji 


( 2 ) 


where the absolute value denotes the undetermined di¬ 
rection of the in-plane DM vector. 

In order to estimate the exchange couplings Ji and J 2 , 
we first analyze magnetic susceptibility based on a weakly 
coupled spin-chain model (J 2 ^ Ji). We reconsider 
the magnetic susceptibility as a function of temperature 
shown in Fig. 3 and fit the data to the result of the high- 
accuracy numerical Bethe ansatz calculations^^ XBethe(^) 
for the ID Heisenberg antiferromagnetic system. A fit of 
the measured magnetic susceptibility (Fig. 5) to 


X(r) =Xo+XBethe(^)’ (3) 

where xo is a temperature-independent susceptibility 
background, yields xo = 1-68(7) x 10“^ cm^/mol Cu, the 
intra-chain coupling Ji = 5.95(2) meV, and the Lande 
^-factor g = 2.16(1). The Curie-Weiss exchange cou¬ 
pling Jew = 12.6 meV obtained from the mean field 
approximation for a number of coordinate spins ^ = 2 is 
substantially larger than Ji obtained from the suscepti¬ 
bility fit to the spin-chain model (| Ji — Jew \ / Jew — 
0.53) suggesting that the coordination number must 
be greater than two. Therefore, the inter-chain in¬ 
teraction is non-negligible. Using Schulz’s calcula¬ 
tion of the inter-chain interaction for parallel chains 
in thle mean field approximation, |Jinter| = 1-^21 = 

Tjv/ (l.28v/ln(5.8Jintra/Tiv)),34 where Antra IS an iutra- 
chain interaction, and applying it to ct-Cu 2 V 207 with 
T]sf = 33.4 K and Antra = Ji = 5.95 meV, we obtain 
IJ 2 I ^ 1.42 meV. Even though the sign of J 2 cannot 
be determined from this calculation, the discrepancy be¬ 
tween Jew for 2 : = 2 and Ji suggests that the inter-chain 
interaction is dominantly antiferromagnetic and hence 
J 2 > 0. Neutron diffraction, which will be discussed 
in the next section, reveals antiferromagnetic alignment 
of spins in the [011] and [Oil] chains on different de¬ 
planes as shown in Fig. 9. Therefore, the leading anti¬ 
ferromagnetic inter-chain interactions (denoted by green 
bonds in Fig. 9(a)), which induce the anti-parallel ar¬ 
rangement of spins along the a-axis, link the chains on 
the adjacent planes giving rise to the helical-honeycomb 
spin network. The Cu-O-Cu bond angle of this inter¬ 
chain interaction, which is greater than 90° (Fig. 1(c)), 
is also consistent with the antiferromangetic coupling. 
The coupling between the parallel chains in the same de¬ 
plane is expected to be much smaller because the bonding 
is not of the Cu-O-Cu type but must be via the V 2 O 7 
double tetrahedra. The value of IJ 2 /A 1 I 0.23 is sub¬ 
stantially larger than that measured in other ID systems 
such as KCuFs, Sr 2 Cu 03 , and BaCu 2 Si 2 07 , where the 
value of I Anter/Aintral IS 0.001 — 0.01,^^“^^ Suggesting that 
the magnetism in a-Cu 2 V 207 is far from being an ideal 
realization of the ID spin system, and probably invalidat¬ 
ing the above analysis as well as other previous studies. 
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TABLE L Parameters obtained from the fit of magnetic susceptibility with H a using different lattice models. 



Bethe ansatz 

calculation^^ 
(coupled spin-chain) 

uniform 

spin-chain 

QMC (this work) 

helical-honeycomb 

Ji > J 2 Ji < J 2 

Ji [meV] 

5.95(2) 

5.95(2) 

5.79(1) 

4.10(1) 

J 2 [meV] 

~1.42 

- 

2.61(1) 

6.31(1) 

9 

2.16(1) 

2.16(1) 

2.24(1) 

2.25(1) 

Xq [10~^ cm^/mol Cu] 

1.68(7) 

1.69(7) 

0.94(6) 

0.89(3) 

\Dp\ [meV] 

0.836(3) 

0.836(3) 

0.814(1) 

0.576(1) 



FIG. 5. (Color online) (a) Fit of magnetic susceptibility with 
iL T a to the Bethe ansatz numerical result (black-dashed 
line), QMC simulations on the uniform spin-chain model 
(green line), QMC simulations on the helical-honeycomb 
model with Ji > J 2 (red), and with Ji < J 2 (blue). (b) Resid¬ 
ual curves (|xdata — XfitI) between the data and fit curves of 
different lattice models. 


which are based on the coupled spin-chain model. 
Furthermore, we note that since the dominating inter¬ 
chain interaction is between the [Oil] and [Oil] chains, 
of which the spin coordination number is different from 
that of the parallel chains, Schulz’s results used above 
could potentially provide an erroneous magnitude of J 2 . 

To obtain the better estimates of Ji and J 2 , we per¬ 
form QMC simulations with the LOOP algorithm on the 
helical-honeycomb spin network, which shows a discon¬ 
nected hexagon when viewed along the 6-axis (Fig. 9(a) 
and (c)). The magnetic susceptibility was computed and 
fit to the data. In addition, the susceptibility for the 


uniform spin-chain model was calculated and compared 
with XBethe(^) (^9- ^)- QMC numerical results were 
first fit using the following formula 

= (4) 

where t = kET/J^a.^ and the Fade approximant is given 
by 

1 + E 

= -. (5) 

1 + y 

n=l 

The coefficient 1/4 is from 5'(5' + l)/3 for S = 1/2 and 
the factor 1/t ensures that x*(t) approaches the Curie 
law at high temperature. Jmax is the maximum of Ji 
and J 2 and the smaller interaction is equal to a Jmax, 
where the parameter a is fixed for each QMC simulation. 
The numerical parameters Nn and Dn were obtained up 
to g = 5 and r = 6 from the fit. The resulting Fade 
approximant for each QMC calculation was then used to 
fit the measured susceptibility data using 

x(^') = Xo+X qmc(^)> (6) 

where 

Xqmc( 2^) = f^^x* (kBT/J^,^), (7) 

to obtain Jmax and g. Na, and ks are the Avogadro 
constant, the Bohr magneton, and the Boltzmann con¬ 
stant, respectively. For the helical-honeycomb model, a 
is manually adjusted until the best fit is acquired. Fig¬ 
ure 5 shows the best fits for both uniform spin-chain 
and helical-honeycomb models where the obtained fit 
parameters are summarized in Table I for comparison. 
The QMC simulation on the spin-chain model exactly 
match the result obtained from the Bethe ansatz calcu¬ 
lations with the same magnitude of Ji, proving the valid¬ 
ity to the QMC simulations. For the helical-honeycomb 
model, two sets of the exchange parameters, for Ji > J 2 
Jmax = 5.79(1) meV, a = 0.45, and g = 2.24(1) and for 
Ji < J 2 Jmax = 6.31(1) meV, a = 0.65, and g = 2.25(1), 
fit the data equally well and provide a noticeably better 
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fit than the spin-chain model for most of the temperature 
range from 300 K down to ^ T/v as shown by the resid¬ 
ual analysis in Fig. 5(b). The obtained average exchange 
interaction J = (2Ji + J 2 ) /3 for the coupled spin-chain 
model is much lower than the mean-field value Jew ob¬ 
tained from the Curie-Weiss fit with | J — Jew \ / Jew — 
0.47. Our proposed helical-honeycomb model shows a 
slightly smaller deviation with |J— Jew \ j Jew — 0.44 
and 0.42 for J\ > J 2 and Ji < J 2 , respectively. Addi¬ 
tional exchange and anisotropic interactions could affect 
the value of J, potentially leading to closer agreement 
between J and Jew- Unfortunately, based on combined 
susceptibility analysis and QMC simulations, we were un¬ 
able to uniquely identify the leading coupling constant. 
It was previously proposed that the leading coupling con¬ 
nects Cul and Cu4 via Ji forming the zigzag chains along 
[Oil] and [Oil] directions (Fig. 1(a)), but the competi¬ 
tion between ferromagnetic coupling via Cul-0(3)-Cu4 
with the bond angle close to 90° and antiferromagnetic 
coupling via Cul-0(2)-Cu4 with the bond angle of 106° 
(Fig. 1(b)) could render a weaker net antiferromagnetic 
interaction for Ji than the non-competing antiferromag¬ 
netic coupling J 2 between Cul and Cu 2 via two equiva¬ 
lent Cul-0(1)-Cu2 bridges (Fig. 1(c)). Further theoret¬ 
ical analyses based on first-principle calculations, which 
will provide complementary support to our analysis, are 
desirable in order to identify the leading interaction. 

Given Ji obtained from the two helical-honeycomb 
models discussed above and the canting angle 77 = 4.0° 
from the magnetization data, we estimate the in-plane 
DM parameter \Dp\ using Eq. (2) to be 0.814(1) meV 
for the helical-honeycomb model with Ji > J 2 and 
0.576(1) meV for Ji < J 2 , or 0.14Ji. The val¬ 
ues of \Dp\ for all considered lattice models are sum¬ 
marized in Table L The magnitude of the DM vec¬ 
tor shows good agreement with Moriya’s calculation^^ 
as /S.g/g ^ 0.12 ^ |Dp/Ji|, where g is the free electron 
Lande ^-factor and denotes its shift caused by the 
crystalline environment. The helical-honeycomb lattice 
formed by Ji and J 2 of comparable strength and the 
presence of the DM interactions in a-Gvi 2 ^ 2^7 induce 
the magnetically ordered state below the Neel tempera¬ 
ture. Hence, the long-range order observed in this system 
does not defy the Mermin-Wagner theorem. 


C. Neutron diffraction 

In order to investigate the magnetic structure of the 
ordered state, we performed elastic neutron scattering 
to search for magnetic Bragg reflections. Figures 6 (b) 
and (c) show clear extra scattering intensity below Tat at 
( 0 , 2 , 0 ), indicative of magnetic scattering; small peaks at 
( 0 , 2 , 0 ) and ( 6 , 0 , 0 ), which are structurally forbidden by 
symmetry, at 50 K are due to higher-order neutron con¬ 
tamination. On the other hand, the intensities at (4,0,0) 
and (6,0,0), which are also allowed magnetic Bragg re¬ 
flections, do not show a significant change below T/v- 



HOO [r.l.u.] 



FIG. 6. (Color online) (a) An intensity map of single-crystal 
q;-Cu 2 V 2 07 was measured in the (/i, /c, 0) plane at 2.5 K. The 
intensity is presented in a log scale, (b) and (c) show the 
0 — 2^-scans around the (4,0,0), (0,2,0), and (6,0,0) Bragg 
reflections at 2.5 K and 50 K, respectively. The solid line is a 
guide to the eye. Error bars represent one standard deviation 
throughout the article. 


In neutron scattering, only the spin component that 
is perpendicular to the momentum transfer contributes 
to the magnetic scattering intensity, due to the dipole- 
dipole interaction as described by the geometric factor 
~ QaQ/s)^^^ where a denotes the spin compo- 

a,l3 

nents and Qa is the unit vector of Q along the compo¬ 
nent a. Eor (a-Gu 2 V 207 most of the spin component is 
parallel to the a-axis, which is also evidenced by the mag¬ 
netic susceptibility. Hence the magnetic Bragg reflections 
(h, 0 , 0 ) become negligibly small. 

The integrated intensity of the (0, 2, 0) magnetic Bragg 
reflection as a function of temperature (Eig. 3 inset) 
shows an upturn around 33 K indicating a transition 
to the antiferromagnetic ordered state, coincident with 
the jump in the magnetic susceptibility (Eig. 3(a)). A 
fit of the order parameter to the power law, I{T) oc 
(1 — T/T/v)^^, for 24 K < T < 34 K gives a critical expo- 





















TABLE IL Refined values of fractional coordinates of a- 
CU 2 V 2 O 7 from single-crystal x-ray diffraction measured at 
room temperature and from single-crystal neutron diffraction 
measured at 50 K. 

Atom x/a y/b zjc U 

x-ray diffraction 

Cu 0.16572(5) 0.3646(1) 0.7545(1) 0.0143(3) 

V 0.19898(5) 0.4055(1) 0.2370(2) 0.0067(3) 

0(1) 0.2455(3) 0.5631(9) 0.274(1) 0.022(1) 

0(2) 0.1445(3) 0.4375(6) 0.0308(8) 0.0100(9) 

0(3) 0.1617(3) 0.3475(7) 0.4575(9) 0.014(1) 

0(4) \ I 0.156(2) 0.030(2) 

Ri = 0.039, wRi = 0.095, GOF = 1.063 
neutron diffraction 

Cu 0.16551(2) 0.36460(5) 0.7520(1) 0.0036(1) 

V 0.1990(3) 0.4046(8) 0.237(2) 0.0006(9) 

0(1) 0.24605(3) 0.56165(8) 0.2723(2) 0.0072(2) 

0(2) 0.14428(3) 0.43776(7) 0.0286(2) 0.0040(2) 

0(3) 0.16200(3) 0.34608(9) 0.4551(2) 0.0052(2) 

0(4) I I 0.1507(3) 0.0087(3) 

Ri = 0.081, wR2 = 0.073, GOF = 1.215 


nent (3 = 0.21(1) and Tn = 33.4(1) K. The fit value of p 
is typical for low-dimensional magnetic systems. 

Figure 7 shows the powder neutron diffraction data col¬ 
lected at 3 K and 40 K at Echidna to investigate intensity 
difference due to magnetic scattering. Figure 7(a) shows 
a small shoulder peak for the 3 K data as indicated by the 
arrow, which corresponds to the (0,2,0) magnetic Bragg 
reflection. In the inset, the intensity difference between 
3 K and 40 K data clearly shows the magnetic Bragg 
reflections at (1,1,1), (0,2,0), and (3,1,1), and negligi¬ 
ble magnetic scattering intensities at (4,0,0) and (6,0,0). 
This absence of magnetic Bragg intensity at the (h, 0, 0) 
reflections is consistent with the single-crystal data dis¬ 
cussed previously. 

The magnetic structure of Q(-Cu 2 V 207 was analyzed by 
the irreducible representation theory. The calculations 
were carried out using the software Baslreps.^^ The de¬ 
composition of the irreducible representations (IRs) for 
Cu^+ ions (165) can be described by 

r = dTi + 3r2 + 3r3 + 3r4. (8) 

The basis vectors for each IR are summarized in Ta¬ 
ble III. The coupled intra-chain Cu^+ spins are between 
Cul-Cu4 along [011], and between Cu2—Cu3 along [011]. 
As discussed above, the antiferromagnetic spin compo¬ 
nent of the Cu^+ ions is predominantly along the a-axis. 
Therefore, among the four possible magnetic models we 
can simply rule out r 2 and Es, which according to the 
Bertaut’s notation^^ give rise to the Ax and Fx config¬ 
urations, respectively; both IRs give ferromagnetic spin 
component along the a-axis. On the other hand, Ti and 
r 4 , which give rise to the Gx and Cx configurations, re¬ 
spectively, result in the antiferromagnetic arrangement 



20 [deg] 



20 [deg] 

FIG. 7. (Color online) (a) Neutron diffraction data of powder 
a-Cu 2 y20r at 3 K show the magnetic Bragg scattering at 
(0, 2, 0) (indicated by an arrow). The inset shows the intensity 
difference with allowed magnetic Bragg reflections indexed, 
(b) The powder neutron diffraction data collected at 3 K are 
refined using FullProf. 


along the a-axis for the spins along the zigzag chains. 
The canted moments, which are in the 5c-plane, of the 
spins in the same chain for Fi and r 4 are parallel. Hence 
the cross product x Sj does not have a component 
along the a-axis, making the DM component along the 
a-axis irrelevant as previously stated. Both Fi and r 4 
can give rise the weak ferromagnetism observed in the 
magnetic susceptibility. However, Fi gives a better fit to 
the powder neutron diffraction data. Figure 7(b) shows 
the full pattern refinements for Fi with the spin compo¬ 
nent along the a-axis rux as the only fit parameter while 
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TABLE III. Magnetic representations and their basis vectors 
for Cul(x, y, z), Cu2(—a:, —y, z), Cu3(x + 1/4, —y + 1/4, z + 
1/4), and Cu4(—a; + 1/4, y + 1/4, z + 1/4) (see Fig. 1). 
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FIG. 8. (Color online) False color representation of the single¬ 
crystal neutron diffraction measured at 4 K on the Lane 
diffractometer Koala. The halo around the center of the image 
is due to scattering from a polycrystalline aluminium crystal 
holder. 


the other two components, i.e. rriy and were fixed 
to zero. The magnetic i?—factors from the refinements 
for Fi and r 4 yield 0.040 and 0.134, respectively, attest¬ 
ing to the validity of Fi over F 4 . The obtained magnetic 
moment along the a-axis rrix is 0.9(2) jj^b for Fi. 

To further confirm the magnetic structure and obtain 
a better estimate of the ordered moment, we performed 
single-crystal neutron diffraction at the Lane diffractome¬ 
ter, Koala. The diffraction data were collected at 4 K and 
50 K to investigate the nuclear and magnetic structure. 
The Lane diffraction pattern measured at 4 K shows dis¬ 
tinct Bragg peaks (Fig. 8 ). The structural parameters for 
the 50 K data were refined against Fdd2 space group us¬ 
ing ShelXle yielding Ri = 0.081 for 998 reflections with 
Fobs < 4cr(Fo6s). The result is in agreement with the 
single-crystal x-ray diffraction refinement (see Table II). 
For the 4 K data, where the system becomes magnet¬ 
ically ordered, the magnetic structure refinements were 
performed using Jana2006^^ The reflections were refined 
against Shubnikov space group Fd'd'2^ which is equiva¬ 




FIG. 9. (Golor online) (a) The magnetic structure of the 
S = 1/2 Gu^^ spins in 0 -GU 2 V 2 O 7 shows the major spin 
component along the a-axis with the antiferromagnetic ar¬ 
rangement. Blue bonds represent Ji while green bonds J 2 . 
(b) The two zigzag chains connected by Ji (J 2 not shown) 
on adjacent 5c-planes are along [Oil] and [011] and are about 
75° with respect to each other. We note that the canted mo¬ 
ments in the 5c-plane are exaggerated for visualization, and 
the drawn spins do not represent the actual canting, neither 
in terms of magnitude nor direction, (c) The connectivity of 
Ji and J 2 gives rise to the helical-honeycomb pattern when 
viewed along the 5-axis. 


lent to the irreducible representation Fi that gives the 
best fit to the magnetic structure in this system. The 
spin components along the b and c axes were fixed to 
zero due to the unresolved spin canting. The ordered 
moment along the a-axis rrix = 0.93(9) /ib was obtained 
with wR = 0.051. This value of the ordered moment, 
which is consistent with the value obtained from powder 
neutron diffraction, is slightly lower than (but close to) 
the expected value of one This discrepancy could 
be a result from the constrained spin component to only 
the a-axis, discarding the spin canting that is not ob¬ 
tainable from the neutron diffraction data. In addition, 
quantum fluctuations might also play a role in reduc¬ 
ing the ordered moment. The antiferromagnetic spin 
structure on the helical-honeycomb lattice is depicted in 
Fig. 9(a). The spins on the parallel zigzag chains fer- 
romagnetically align as a result of the antiferromagnetic 
coupling J 2 around the helical-honeycomb loop as shown 
in Fig. 9(c), allowing the ordering to propagate along the 
transverse directions and prompting the 3T) long-range 
order. Due to the very small spin canting we cannot de¬ 
termine the magnitude and direction of the canted mo¬ 
ments in the 6 c-plane from the neutron diffraction data. 
However, Fi allows the canted moments of the spins on 
the same chain to be parallel. For the spins on the differ- 
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ent chains that are located on the adjacent 6 c-plane, the 
c-component (^-component) is parallel (anti-parallel) as 
shown in Fig. 9(b). 


IV. SUMMARY 

We have proposed a new spin model to describe mag¬ 
netic properties of (a-Cu 2 V 207 . Combined studies of 
magnetization, QMC simulations, and neutron diffrac¬ 
tion show the helical-honeycomb pattern of spin network 
connected by two different exchange couplings. Magnetic 
susceptibility shows a broad peak at ^ 50 K, which is 
an evidence of rising short-range spin correlations, fol¬ 
lowed by an abrupt increase indicative of a phase transi¬ 
tion to a magnetically ordered state at T/v = 33.4(1) K. 
The Bethe ansatz calculations for the S = 1/2 uniform 
Heisenberg chain fit the H ± a data very well above T/v 
but our proposed helical-honeycomb model shows sub¬ 
stantial improvement of the fit for Ji : J 2 = 1 : 0.45 
with Ji = 5.79(1) meV and for Ti : T 2 = 0.65 : 1 
with J 2 = 6.31(1) meV. Therefore, we conclude that the 
helical-honeycomb model describes the underlying spin 
network of a-GM 2 \f 2^7 more accurately than the previ¬ 
ously held spin-chain model. 

The anisotropy below Tn suggests that the majority 
of the spin component is along the crystallographic a- 
axis, which is confirmed by neutron diffraction exper¬ 
iments. The weak ferromagnetism is a result of spin 
canting within the 6 c-plane due to the DM interactions. 
Magnetization measurements with H _L a show the spon¬ 
taneous magnetization from which the canting angle r] of 
4.0° and the in-plane DM parameter \Dp\ 0.14Ji are 
obtained. 

The analysis of the neutron diffraction data shows that 
the S = 1/2 Cu^+ spins antiferromagnetically align along 
the helical-honeycomb loops with the ordered moment of 
0.93(9) /iB predominantly along the crystallographic a- 
axis. The spin network of two comparable exchange cou¬ 
plings forming the helical-honeycomb lattice and the DM 
interactions lead to the long-range magnetic ordering be¬ 


low T/v- However, due to the complex exchange pathways 
and the presence of weak frustration, the exact value of 
both exchange interactions could deviate from our anal¬ 
ysis. Further theoretical analyses based on first-principle 
calculations and studies of spin dynamics by means of in¬ 
elastic neutron scattering are therefore required in order 
to better determine the exchange interactions and con¬ 
firm the helical-honeycomb spin model. With the current 
availability of large single crystals, an inelastic neutron 
scattering study is possible and will reveal a complete 
picture of the relevant microscopic Hamiltonian param¬ 
eters, as well as influences of the low spin-coordination 
number and quantum fluctuations on spin dynamics. 

Note added in proof After the submission of this 
manuscript, we became aware of the work by Sannigrahi 
et in which they performed density functional the¬ 
ory to calculate the relevant exchange interactions in a- 
CU 2 V 2 O 7 . However, the QMC calculations of the sus¬ 
ceptibility show that their lattice model with the dom¬ 
inant third-nearest-neighbor antiferromagnetic interac¬ 
tions (J 3 = 13.61meV) gives a broad peak at a higher 
temperature than ^ 50 K, inconsistent with the experi¬ 
mental data. 
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